%% FigS23A
% Change i = 9,10,11 for the figures in FigS23A
for i = 11
    
T = i;

load(strcat('C:\Users\AlexCui\Dropbox (Harvard University)\Kim_Group_Harvard\Twist Junctions\2022-05-15 - Interface 66\Cseries2\C00',num2str(T),'Icpos.mat'));

load(strcat('C:\Users\AlexCui\Dropbox (Harvard University)\Kim_Group_Harvard\Twist Junctions\2022-05-15 - Interface 66\Cseries2\C00',num2str(T),'Icneg.mat'));

figure(233)
Ic1 = histogram(Ic5Kpos./0.0018,[400:2:700]); hold on;
Ic2 = histogram(-Ic5Kneg./0.0018,[400:2:700]);hold off



end

xlim([400 700]);
ylim([0 2500]);
hold off

%% FigS23B
% Change i = 1,2,3 for the figures in S23B
for i = 3
    
T = i;

load(strcat('C:\Users\AlexCui\Dropbox (Harvard University)\Kim_Group_Harvard\Twist Junctions\2022-09-17 Interface 74\Fastsweeps\Oseries\O00',num2str(T),'bIcpos.mat'));

load(strcat('C:\Users\AlexCui\Dropbox (Harvard University)\Kim_Group_Harvard\Twist Junctions\2022-09-17 Interface 74\Fastsweeps\Oseries\O00',num2str(T),'bIcneg.mat'));

figure(233)
Ic1 = histogram(Ic5Kpos./0.0018,[700:4:1300]); hold on;
Ic2 = histogram(-Ic5Kneg./0.0018,[700:4:1300]);hold off

end

xlim([700 1300]);
ylim([0 400]);
hold off